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Abstract 

The  objective  of  this  work  is  to  investigate  the  effect  of  clamping  force  on  the  interfacial  contact  resistance  and  the  porosity  of  the  gas  diffusion 
layer  (GDL)  in  a  proton  exchange  membrane  fuel  cell  (PEMFC).  An  optimal  rib  shape  for  the  bipolar  plate  is  developed  to  analyze  the  electrical 
contact  resistance.  We  found  that  the  electrical  contact  resistance  is  determined  by  both  the  clamping  force  and  the  contact  pressure  distribution.  A 
minimum  contact  resistance  can  be  obtained  in  the  case  of  a  constant  contact  pressure  distribution.  The  porosity  of  the  GDLs  underneath  the  rib  of 
the  bipolar  plate  decreases  with  increasing  the  clamping  force,  and  the  void  volume  is  changed  with  the  deformation  of  the  GDLs.  It  is  found  that 
there  exists  an  optimal  rib  width  of  the  bipolar  plates  to  obtain  a  reasonable  combination  of  low  interfacial  contact  resistance  and  good  porosity 
for  the  GDL. 

©  2006  Published  by  Elsevier  B.Y. 
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1.  Introduction 

The  basic  working  mechanism  of  electrochemical  fuel  cells 
(FC)  is  to  convert  fuel  and  oxidant  directly  into  electricity  energy. 
Since  1839,  when  Grove  [1]  demonstrated  a  primitive  fuel  cell, 
great  progress  has  been  made  in  the  research  and  development 
of  fuel  cells.  Today,  fuel  cells  are  considered  to  have  a  great 
future  due  to  their  zero  pollution,  high  efficiency  and  very  low 
noise.  Among  the  fuel  cells  under  development,  the  polymer 
electrolyte  membrane  fuel  cell  (PEMFC)  is  considered  as  the 
first  candidate  for  a  power  source  for  automobiles  because  of 
its  low-temperature  operation,  high  power  density,  fast  start-up 
and  system  robustness.  During  the  past  20  years,  much  attention 
has  been  placed  on  the  research  and  development  of  the  PEMFC 
with  a  Nafion®  membrane  as  electrolyte.  To  increase  the  power 
density,  many  researchers  have  focussed  their  attention  on  novel 
membrane  electrolytes,  catalysts  and  structures,  as  well  as  on 
electrode  materials  (gas  diffusion  layers)  and  their  preparation. 
In  a  single  PEMFC,  a  polymer  membrane  is  inserted  between 
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two  gas  diffusion  layers  (GDLs)  coated  with  catalyst,  which 
is  further  sandwiched  between  two  bipolar  plates  housing  the 
flow  channels.  In  a  practical  application,  multiple  fuel  cells  are 
stacked  together  to  give  a  sufficiently  high  power  and  desired 
voltage. 

The  gas  diffusion  layer,  catalyst  and  membrane  are  usu¬ 
ally  united  together  as  a  membrane  electrode  assembly  (ME A) 
before  stack  fabrication.  During  the  clamping  process  of  a  fuel 
cell  stack,  ME  As,  bipolar  plates,  end  plates  and  gaskets  are 
finally  clamped  together,  as  shown  in  Fig.  1.  The  gas  diffu¬ 
sion  layers  are  made  of  highly  porous  materials,  and  thus  are 
deformed  very  much  in  the  clamping  process.  The  membrane 
may  be  also  deformed  greatly  due  to  its  low  elastic  modulus.  The 
compression  and  deformation  of  the  ME  A  depends  on  many 
parameters,  such  as  the  elastic  modulus  and  thickness  of  the 
MEA,  clamping  force,  the  pressures  of  the  fuel  and  oxidant, 
gasket  material  and  thickness,  temperature,  humidity,  etc.  The 
deformation  of  those  materials  affect  directly  the  efficiency  of 
the  fuel  cell  stack  [2,3].  If  the  clamping  force  is  unreasonably 
small,  the  interfacial  contact  resistance  between  the  gas  diffu¬ 
sion  layer  and  bipolar  plate  will  be  increased  and  thus  the  system 
efficiency  is  reduced  [3-5].  However,  if  the  clamping  force  is 
unreasonably  large,  the  efficiency  will  be  also  decreased.  In 
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addition,  if  it  is  too  large,  the  MEA  components  may  be  dam¬ 
aged,  causing  cell  leakage  and  internal  shorts.  Experiments  have 
showed  that  there  is  an  optimum  clamping  force  for  a  real  fuel 
cell  stack  [3]. 

After  being  clamped  together  in  a  fuel  cell  system,  the  gas 
diffusion  layer  may  be  one  of  the  components  with  the  most 
deformation.  Thus,  the  permeability  or  the  pore  size  may  be 
changed  very  much.  This,  in  turn,  influences  greatly  the  power 
density  or  efficiency  of  the  fuel  cell  system  [6-8],  especially  in 
an  inter-digital  flow  channel  [9]  or  with  serpentine  flow  channels 
[8] .  In  the  present  paper,  we  first  theoretically  study  the  influence 
of  the  clamping  force  and  the  width  of  the  rib  on  the  interfacial 
contact  resistance  and  porosity  of  the  gas  diffusion  layer,  then, 
we  will  discuss  their  possible  effects  on  performance  of  a  fuel 
cell. 

2.  Contact  resistance  model 

Multiple  fuel  cells  are  stacked  together  in  most  practical 
applications  to  provide  sufficiently  high  power  and  voltage.  In 
PEMFCs  a  gas  diffusion  layer  contacts  directly  with  the  bipo¬ 
lar  plates.  The  interfacial  electrical  contact  resistance  between 
them  is  one  of  the  important  factors  leading  to  operational  volt¬ 
age  loss,  and  can  significantly  affect  the  performance  of  a  fuel 
cell  system  [10] .  The  contact  resistance  of  a  pair  of  solid  surfaces 
is  governed  by  the  multi-scale  surface  topography  (roughness) 
and  the  surface  chemical  property  (a  thin  layer  of  oxide  or  con¬ 
tamination).  When  two  solid  surfaces  contact  each  other  with 
a  very  small  compression  force,  only  the  tips  of  the  roughness 
come  into  real  contact  and  thus  the  real  contact  area  is  very 
small  [1 1,12].  It  was  found  that  the  contact  resistance  decreased 
non-linearly  with  compression  force. 

Several  groups  of  experimental  data  on  the  interfacial  contact 
resistance  of  a  single  fuel  cell  or  the  total  resistance  that  includes 
the  bulk  resistance  are  available  in  the  literature  [5,13,14].  In 
a  fuel  cell,  the  contact  resistance  may  totally  differ  from  that 
of  elastic  solids  since  the  GDL  is  a  layer  of  composite-porous 
material.  Mishra  et  al.  [13]  measured  the  total  resistance  and 
the  interfacial  contact  resistance  between  different  gas  diffu¬ 
sion  layers  and  graphite  bipolar  plates.  A  fractal  asperity-based 
model  was  adopted  to  predict  the  interfacial  contact  resistance 
as  a  function  of  pressure,  material  properties  and  surface  geom¬ 
etry.  Wang  et  al.  [5]  reported  the  interfacial  contact  resistances 
between  a  Toray  carbon  paper  (Electro-Chem  Inc.)  and  four 


types  of  stainless  steel  bipolar  plates.  A  stainless  steel  bipolar 
plate  is  considered  to  be  one  of  the  best  candidates  for  bipo¬ 
lar  plates  due  to  its  low  cost,  high  strength,  ease  of  machining 
and  shaping  into  thin  sheets,  as  well  as  its  high  corrosion  resis¬ 
tance  [5].  Ihonen  et  al.  [14]  conducted  in  situ  measurements  of 
the  interfacial  contact  resistance  as  a  function  of  time,  clamping 
pressure,  gas  pressure  and  current  density.  All  the  experimental 
studies  mentioned  above  are  focused  on  the  effects  of  clamp¬ 
ing  pressure  on  the  interfacial  contact  resistance.  It  was  found 
that  the  interfacial  contact  resistance  decreases  quickly  with  the 
clamping  force  at  a  small  clamping  force,  and  reaches  almost  a 
saturation  level  at  high  clamping  force. 

There  are  many  researches  on  the  electricity /heat  conduction 
of  two  contact  surfaces  [11,12,15],  giving  a  few  topography- 
based  theoretical  models.  The  mechanisms  of  electricity  and 
heat  conduction  are  similar  to  each  other,  strongly  depending 
on  the  real  contact  area.  However,  the  topography  parameters 
used  in  the  model  are  very  complex  and  difficult  to  measure.  In 
addition,  no  research  work  was  found  to  consider  such  a  special 
structure  of  random  fibre  network  as  a  carbon  paper.  Different 
fabrications  of  the  GDLs  and  bipolar  plates  produce  different 
interfacial  contact  resistances.  In  the  present  work,  we  use  an 
equation  based  on  the  work  of  Mishra  et  al.  [13]  to  calculate 
the  total  interfacial  contact  resistance.  Generally  speaking,  the 
relationship  between  interfacial  resistance  and  contact  pressure 
can  be  expressed  by  following  equation: 


where  Aa  is  the  apparent  contact  area  of  the  interface,  p  the 
contact  pressure  and  A,  B  and  C  are  parameters  determined  by 
experiments.  According  to  Eq.  (1),  we  can  get  the  interfacial 
contact  resistivity: 


We  found  that  the  interfacial  contact  resistivity  measured  by 
Mishra  et  al.  [13]  and  Wang  et  al.  [5]  can  be  fitted  by  Eq.  (2) 
with  a  pretty  good  accuracy,  as  shown  in  Fig.  2.  The  values  of 
fitted  parameters  are  shown  in  Table  1 .  A  high  contact  pressure 
can  decrease  the  total  electrical  contact  resistance.  On  the  other 
hand,  the  clamping  force  can  change  the  porosity  of  gas  diffusion 
layer.  A  higher  porosity  may  be  expected  to  decrease  the  mass 
transfer  resistance  of  a  gas  phase,  especially  in  a  PEM  fuel  cell 
with  inter-digital  flow  field  as  shown  in  Fig.  3. 

Gas  diffusion  layer  does  much  more  than  diffusing  the  gas 
and  forming  an  electrical  connection  between  the  catalyst  layer 
and  the  bipolar  plate.  In  addition,  it  carries  water,  an  important 
reaction  product,  away  from  the  electrolyte  surface,  and  also 
forms  a  protective  layer  over  the  very  thin  layer  of  the  catalyst. 
When  a  fuel  cell  system  is  stacked,  the  void  volume  in  gas  dif¬ 
fusion  layer  underneath  the  rib  will  decrease  with  the  clamping 
force.  For  each  fuel  cell  system,  there  might  be  an  optimal  choice 
of  the  clamping  force. 
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Table  1 

Parameters  of  the  interfacial  contact  resistance  function  for  gas  diffusion  layers  and  the  bipolar  plate 
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Bipolar  plate 

GDL 

A  (m<2  cm2) 

B  (MPa) 

C 

Reference 

Graphite 

GDL- 1  OB  A 

3.32 

1.01 

0.534 

Mishra  et  al.  [13] 

Graphite 

GDL-10BB 

3.72 

0.966 

0.692 

Mishra  et  al.  [13] 

Stainless  steel 

Carbon  paper 

81.4 

2.52 

1.07 

Wang  et  al.  [5] 

GDL-10BA  and  GDL-10BB  are  two  kinds  of  paper-based  GDL  (SGL  Carbon  Group). 


Fig.  2.  Variation  of  the  interfacial  contact  resistance  over  a  range  of  clamping 
pressures. 


Gas  out 


2 

3 


X - \ 


(b) 

Fig.  3.  (a)  Schematic  of  an  inter-digital  flow  field  and  (b)  a  typical  cross-section. 
1,  Bipolar  plate  (current  collector);  2,  inlet  gas  channel;  3,  outlet  gas  channel; 
4,  gas  diffusion  layer;  5,  catalyst  layer  and  6,  membrane. 


3.1.  Finite  element  meshes  and  material  parameters 

We  used  the  finite  element  method  (FEM)  to  analyze  the  con¬ 
tact  pressure  and  the  bulk  strain  of  the  gas  diffusion  layer.  The 
FEM  analysis  model  is  shown  in  Fig.  4.  In  order  to  save  compu¬ 
tation  time,  the  computational  domain  consists  of  half  a  period 
of  the  flow  channel  and  rib,  i.e.,  the  left  half  of  the  area  of  Fig.  3b. 
It  is  assumed  that  the  structures  of  the  cathode  and  anode  of  the 
bipolar  plates  have  the  same  shape,  and  the  compressive  modu¬ 
lus  of  catalyst  layer  and  membrane  are  much  higher  than  that  of 
the  gas  diffusion  layer.  The  round  radius  of  the  rib  will  change 
the  local  pressure  distribution,  and  it  is  very  important  to  avoid 
stress  concentration  leading  to  possible  damage  of  the  local  strc- 
tural  materials.  The  normal  displacements  at  the  left,  right  and 
bottom  boundaries  of  the  gas  diffusion  layer  were  restricted  due 
to  the  symmetrity  of  the  structure.  The  vertical  displacement  of 
the  rib  is  given  downwards,  simulating  the  stacking  process. 

The  problem  mentioned  above  can  be  described  as  a  large 
deformation  contact  problem  with  plane-strain  and  free  bound¬ 
ary.  The  meshes  of  finite  elements  were  designed  as  shown  in 
Fig.  4.  Plane-strain  elements  with  eight  nodes  were  used  and 
contact  elements  were  built  in  the  contact  region.  Table  2  shows 


Table  2 


Compression  modulus  as  a  function  of  pressure  for  GDLs[13] 


GDL 

Pressure  range  (MPa) 

E  (MPa) 

GDL- 1 0B  A 

0.00-0.15 

1.71 

0.15-1.12 

4.59 

0.12-3.00 

6.16 

GDL-10BB 

0.00-0.16 

1.81 

0.16-0.52 

5.11 

0.52-3.00 

8.57 
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the  compressive  modulus  of  the  gas  diffusion  layers  as  a  func¬ 
tion  of  the  relevant  pressure  range  [13].  It  was  observed  that 
the  GDLs  have  three  regions  of  linear  elasticity.  The  GDL  is  a 
layer  of  anisotropic  porous  material.  The  elastic  modulus  in  v 
direction  (in-plane  direction)  is  different  from  that  in  y  direc¬ 
tion  (normal  direction).  However,  due  to  the  symmetry  of  the 
fuel  cell  structure  and  the  very  thin  GDL  layer  compared  to 
the  other  dimensions,  the  in-plane  deformation  was  constrained 
in  the  numerical  analysis.  In  fact  the  in-plane  deformation  has 
relatively  insignificant  effects  on  the  contact  pressure  and  bulk 
strain.  Consequently,  the  in-plane  mechanical  properties  are 
not  important  to  the  parameters  studied  in  the  present  paper, 
and  we  assume  that  the  GDL  is  isotropic.  In  order  to  check 
the  validity  of  the  assumption,  we  used  the  in-plane  elastic 
modulus  ranging  from  5  to  50  MPa,  and  found  that  this  gave 
rise  to  a  change  of  the  contact  pressure  and  bulk  strain  less 
than  1.5%. 


3.2.  Interfacial  contact  resistance 

The  interfacial  contact  resistance  between  two  surfaces  is 
considered  as  an  electric  circuit  in  a  parallel  connection,  as 
shown  in  Fig.  5,  supposing  that  the  surface  is  equipotential.  The 
interfacial  contact  resistance  can  be  calculated  using  the  equa¬ 
tion  given  in  Fig.  5  with  the  continuous  distribution  of  interfacial 
contact  resistance,  where  p  is  the  electrical  resistivity,  R b  the 
bulk  resistance  of  bipolar  plate,  Rq  the  bulk  resistance  of  gas 
diffusion  layer  and  7?g/b  is  the  contact  resistance  between  gas 
diffusion  and  bipolar  plate.  After  the  numerical  solution  of  the 
interfacial  contact  pressure,  we  can  get  the  discrete  contact  pres¬ 
sure  distribution.  The  total  interfacial  contact  resistance  can  be 
calculated  using  the  following  equation: 

R  =  X  =  5/  ^ 

2^i=i  Rt  2^i= l  A(B/pinf 


where  Si  and  n  denote  the  boundary  area  of  the  ith  contact  ele¬ 
ment  and  the  number  of  contact  elements,  respectively.  The  total 
compression  force,  F,  has  the  following  relationship  with  the 
contact  stress  components: 

n 

F  =  +  p\x  (4) 

i=  1 

where  pln  is  the  normal  contact  stress  component  (contact  pres¬ 
sure),  p\  the  tangential  contact  stress  component,  n  and  r  are  the 
unit  vectors  in  normal  and  tangential  directions,  respectively. 


3.3.  Porosity  of  the  GDL 


The  GDL  is  a  porous  elastic  material.  The  porosity  of  porous 
media  is  defined  as: 


(5) 


where  Vp  is  the  void  volume  and  V  is  the  total  volume  of  porous 
media.  Based  on  the  definition  of  porosity,  we  can  get  the  poros¬ 
ity  of  GDL  after  stacking: 


Vp  Vo<po  -  Vo£v  <P0  -  £v 

(D  =  -=?-  =  -  =  -  (6) 

V  Vo(l  -  ev)  1  -  ev 

where  Vo  is  the  initial  total  volume  of  porous  media,  (po  the  initial 
porosity  and  £v  is  the  bulk  strain  (ev  =  £x  +  £y  +  sz).  The  top-bar 
means  the  variables  after  stacking.  We  also  defined  a  relative 
porosity  for  comparison  with  the  initial  void  volume: 


Vo^o  -  V0£v 
Vo 


—  (p0  £v 


(7) 


Eqs.  (6)  and  (7)  are  based  on  the  small  deformation  assump¬ 
tion,  and  the  error  is  about  5.7%  when  the  bulk  strain  is  within 
0.3.  Obviously,  the  porosity  of  a  gas  diffusion  layer  after  stack¬ 
ing  loading  is  not  uniform  any  more.  The  minimum  porosity  of 


T 


Bipolar  plate 
Gas  diffusion  layer 


rG/B 


Fig.  5.  Schematic  of  the  contact  electrical  circuit  in  a  parallel  connection. 
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gas  diffusion  layer  occurs  underneath  the  middle  of  rib  (x  =  0  in 
Fig.  4),  where  the  GDL  gives  a  largest  deformation. 

4.  Numerical  results  and  analysis 

4.1.  Optimization  of  the  contact  pressure  distribution 

From  Eq.  (3),  the  total  interfacial  contact  resistance  for  dif¬ 
ferent  contact  pressure  distributions  may  be  totally  different, 
even  with  the  same  clamping  force.  Distributions  of  the  contact 
pressures  are  determined  by  the  mechanical  and  geometrical 
properties  of  the  bipolar  plates.  An  optimum  rib  shape  produc¬ 
ing  a  minimum  interfacial  contact  resistance  is  expected  in  fuel 
cell  stack  design.  However,  it  is  difficult  to  exactly  optimize 
the  rib  shape  directly.  The  first  question  is  what  a  contact  pres¬ 
sure  distribution  gives  a  minimum  contact  resistance  when  the 
total  compression  force  is  fixed.  The  optimum  problem  can  be 
described  as  follows: 

find  pn(x) 

min  .  R  =  R(pn ) 

/  W 

s.t.  /  pndx  =  F 

Jo 


The  design  variable  is  the  contact  pressure  distribution,  and 
the  objective  function,  i.e.,  contact  resistance,  is  a  function  of 
contact  pressure  distribution.  In  the  present  paper,  the  optimum 
contact  pressure  was  first  obtained,  and  then  the  optimum  rib 
shape,  which  gives  an  optimum  contact  pressure  distribution, 
was  chosen.  The  optimization  process  shows  that  a  constant 
pressure  distribution  is  an  optimum  pressure  distribution: 

P(x)  =  ^  (9) 

where  L  is  the  length  of  contact  region.  The  second  question 
is  what  a  rib  shape  of  the  bipolar  plate  gives  a  constant  pres¬ 
sure  distribution.  The  pressure  distribution  is  almost  a  constant 
value  when  the  cross-shape  of  rib  is  rectangular.  When  the  GDL 
is  loaded  with  a  constant  distributed  pressure,  an  optimum  rib 
shape  can  be  found  approximately  from  the  deformation  of  the 
GDL  surface  (see  Fig.  6a),  similar  to  a  rectangular  rib  with  round 
comer.  However,  the  optimum  rib  may  be  more  difficult  for 
machining  than  a  traditional  rectangular  rib. 

The  contact  pressure  approaches  to  infinity  if  the  corner  of 
rectangle  rib  is  a  right  angle,  and  an  increase  of  the  round  radius 
reduces  the  local  pressure.  We  should  design  a  reasonable  round 
radius  to  reduce  the  local  stress  concentration.  However,  too  big 


L=0.50mm 


L=0.75mm 


L=1.00mm 


(a) 


(c)  Width  of  rib,  mm 


Fig.  6.  (a)  Optimum  rib  shape  with  different  widths  of  rib  when  A.  =  2  mm,  hg  =  0.5  mm  and  the  contact  pressure  is  1  MPa;  (b)  interfacial  contact  resistance  at  various 
rib  widths  and  clamping  forces  for  graphite  bipolar  plate  and  (c)  interfacial  contact  resistance  at  various  rib  widths  and  clamping  forces  for  stainless  steel  bipolar 
plate. 
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a  round  radius  will  increase  the  contact  resistance.  The  reason 
is  shown  later  from  the  comparison  of  the  contact  resistances 
produced  by  a  rectangular  rib  and  a  semicircular  rib.  In  this 
study,  the  round  radius  is  0.04  mm. 

4.2.  Effects  of  rib  length  and  clamping  force  on  the 
interfacial  contact  resistance 


Using  the  experiential  formula  (3),  one  obtains: 


=  A 


c 


L 


c-i 


From  Eq.  (10),  one  obtains: 


dR 

dL 


lc~2(c  —  i) 


(10) 


(ii) 


The  experimental  observations  [13]  show  that  C  is  a  constant 
less  than  1  for  a  graphite  bipolar  plate.  Therefore,  the  interfa¬ 
cial  contact  resistance  decreases  with  increasing  the  rib  width. 
However,  it  is  greater  than  1  for  stainless  steel.  The  interfa¬ 
cial  contact  resistance  was  calculated  at  various  rib  widths,  as 
shown  in  Fig.  6b.  It  is  observed  that  the  interfacial  contact  resis¬ 
tance  decreases  rapidly  at  a  small  rib  width  and  then  decreases 
gradually  at  a  large  rib  width,  and  the  resistance  of  GDL-10BA 
with  lower  hardness  decreases  more  quickly  than  GDL-10BB. 
When  the  rib  width  is  fixed,  the  interfacial  contact  resistance 
is  sensitive  to  the  clamping  force,  especially  at  a  value  smaller 
than  0.6  kN  m-1 .  For  a  stainless  bipolar  plate,  C  is  greater  than 
1,  as  a  result  the  value  of  Eq.  (11)  is  positive,  and  the  interfa¬ 
cial  contact  resistance  increases  with  increasing  rib  width  (see 
Fig.  6c).  They  are  52.1  and  60.8  m^m  at  the  rib  width  of  0.2 
and  1.8  mm  (k  =  2  mm),  respectively,  when  the  clamping  force 
is  0.4  kNm-1. 


(a)  Clamping  force,  kN/m 


Fig.  7.  Porosity  (a)  and  relative  porosity  (b)  decrease  as  a  function  of  clamping 
force  and  rib  width. 


4.3.  Porosity  distribution  in  the  GDL 

Now  we  will  discuss  the  effect  of  clamping  force  on  the  inter¬ 
facial  contact  resistance.  On  the  one  hand,  a  higher  clamping 
force  and  larger  rib  width  are  helpful  to  reduce  the  total  interfa¬ 
cial  contact  resistance  for  a  graphite  bipolar  plate.  On  the  other 
hand,  a  higher  clamping  force  results  in  a  decrease  of  porosity 
(i.e.,  the  void  volume  in  gas  diffusion  layer)  and  an  increase  of 
the  gas  transport  resistance.  At  the  same  time,  electrode  flood¬ 
ing  may  occur.  The  porosity  is  a  function  of  the  rib  width  and 
clamping  force  (Fig.  7).  It  is  observed  that  the  relative  porosity, 
defined  in  Eq.  (7),  decreases  much  faster  than  the  porosity  with  a 
clamping  force.  This  indicates  that  the  deformation  of  the  GDL 
and  clamping  force  have  to  be  considered  when  we  analyze  the 
gas  transport. 

Most  of  the  studies  in  the  development  of  a  theoretical  model 
for  PEM  fuel  cells  have  used  for  simplicity  the  assumption  of 
constant  porosity  in  the  GDL  [16-18].  But  Roshandel  et  al.  [19] 
assumed  a  variation  of  porosity  across  the  GDL  as: 


where  £°  is  the  initial  porosity  and  £comp  is  the  porosity  of  the 
GDL  after  stacking.  These  models,  however,  cannot  reflect  the 
importance  of  the  void  volume  change  caused  by  a  clamping 
force. 

4.4.  Comparison  of  different  GDLs  and  rib  shapes 

Two  kinds  of  gas  diffusion  layers  (GDL- 1 0B A  and  GDL- 
10BB)  with  different  initial  porosity  and  compressive  modulus 
are  studied.  Effects  of  clamping  force  on  contact  resistance  and 
relative  porosity  are  shown  in  Fig.  8  (the  width  of  rib  is  1  mm).  It 
is  observed  that  there  is  an  intersection  of  the  two  relative  poros¬ 
ity  curves.  At  the  lower  clamping  force  area,  the  GDL- 1 0B A 
with  smaller  contact  resistance  and  larger  relative  porosity  gives 
rise  to  a  better  performance  than  GDL-10BB.  For  the  higher 
clamping  force  area,  however,  the  contact  resistances  are  almost 
the  same  for  both  of  the  GDLs,  but  the  relative  porosity  of  the 
GDL-10BB  is  larger  than  that  of  the  GDL-10BA.  Thus,  the 
GDL-10BB  may  give  rise  to  a  better  performance  than  GDL- 
10BA. 

Based  on  the  foregoing  discussion,  we  find  that  the  clamping 
force  affects  the  interfacial  contact  resistance  and  the  porosity 
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Clamping  force,  kN/m 


Fig.  8.  Effects  of  clamping  force  on  contact  resistance  and  relative  porosity. 


Fig.  10.  Effect  of  rib  shape  on  contact  resistance  and  porosity. 


in  two  opposite  ways.  A  high  clamping  force  can  give  a  low 
interfacial  contact  resistance  as  expected,  but  also  gives  rise  to  a 
low  porosity  of  GDL,  which  may  decrease  the  system  efficiency. 
The  shape  of  the  rib  is  one  of  the  important  factors  affecting  the 
contact  pressure  distribution,  contact  resistance  and  the  porosity 
of  GDL.  We  now  study  the  effects  of  rib  shape  (rectangular  and 
semicircular)  on  the  contact  resistance  and  porosity  of  GDLs. 
Fig.  9  shows  distribution  of  the  porosity  with  the  same  clamping 
force  but  different  shapes  of  rib  cross-sections.  It  is  observed 
that  the  porosity  of  the  GDL  contacting  with  the  semicircular 
rib  is  smaller  than  that  when  contacting  with  a  rectangular  rib. 


The  permeability,  which  is  related  to  the  porosity,  controls  the 
fuel  flow  distribution  in  fuel  cells.  The  fuel  flow  distribution  has 
an  important  influence  on  the  power  density  [20-22].  As  can 
be  seen  from  Fig.  10,  the  rib  with  a  semicircular  cross-section 
results  in  a  higher  contact  resistance  and  lower  porosity  than  the 
rib  with  rectangle  cross-section.  When  the  clamping  force  per 
unit  length  ranges  from  0.2  to  ~  1 .0  kN  m-1 ,  the  rectangular  rib 
is  better  than  the  semicircular  rib. 

5.  Discussions 


GDL-10BA 
>i=2  mm 
L=1mm 
F=750N/m 
00=0.88 


0.765  0.779  0.792  0.806  0.820  0.834  0.847  0.861  0.875  0.889 
(a) 


GDL-10BA 
7=3mm 
L=1  mm 
F=750N/m 
r=0.04mm 
00=0.88 


0.789  0.800  0.811  0.821  0.832  0.843  0.854  0.865  0.876  0.886 

(b) 

Fig.  9.  The  porosity  distribution  in  the  GDL  for  different  cross-sections  of  rib: 
(a)  semicircular  and  (b)  rectangular  with  a  round  corner. 


Based  on  the  analysis  mentioned  above,  it  is  observed  that 
the  geometric  parameters  of  bipolar  plate,  the  rib  shape,  the 
physical  properties  of  GDL  and  bipolar  plate,  and  the  clamping 
force  will  influence  the  electric  resistance  and  the  porosity  in 
GDLs.  Experiments  of  Hentall  et  al.  [23]  and  Lee  et  al.  [10] 
indicated  that  changes  in  electric  resistance  and  the  porosity  in 
the  GDL  affect  the  power  density  of  a  fuel  cell  system.  In  fact, 
there  are  more  complex  parameters,  such  as  the  fuel  pressure 
inside  the  fuel  cell,  temperature,  humidity,  thickness  of  gasket, 
working  together  to  affect  the  resistance  and  porosity.  What’s 
more,  all  the  parameters  may  influence  to  each  other. 

In  the  present  paper,  the  compression  force  between  a  GDL 
and  a  bipolar  plate  is  found  to  be  one  of  the  key  factors  to 
affect  the  electric  resistance  and  porosity  of  the  GDL.  We  should 
notice  that  in  a  real  PEMFC  system  the  fuel  pressure  supplied, 
mechanical  and  geometrical  properties  of  the  gaskets  and  the 
contact  pressure  between  GDLs  and  bipolar  plates  work  together 
to  determine  the  total  clamping  force,  as  shown  in  Fig.  11. 
There  are  several  questions  which  should  be  paid  more  atten¬ 
tion.  The  volume  of  the  channel  is  changed  more  than  10% 
due  to  the  deformation  of  the  GDL  when  the  nominal  clamping 
pressure  (total  clamping  force  divided  by  the  area  of  GDL)  is 
about  0.4  MPa.  The  effect  of  change  in  the  channel  volume  on 
the  power  density  of  the  fuel  cell  system  should  be  considered 
according  to  the  work  of  Kumar  and  Reddy  [24].  Therefore, 
the  effect  of  the  clamping  force  on  the  power  density  may  be 
deduced  from  the  change  in  channel  volume  besides  the  electric 
resistance  and  porosity.  On  the  other  hand,  the  total  clamping 
force  is  affected  by  many  parameters,  such  as  the  inner  gas  pres- 
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-  Bipolar  plate 


Gasket 


layer  will  decrease  with  increasing  clamping  force,  especially 
with  a  semicircular  rib.  For  a  stainless  bipolar  plate  there  seems 
to  be  an  optimal  rib  width  to  give  a  reasonable  low  contact  resis¬ 
tance  and  a  high  porosity. 

The  constant  pressure  distribution  is  an  optimum  distribu¬ 
tion  for  getting  minimum  interfacial  electric  resistance.  This  is 
a  preparatory  work  to  optimize  the  fuel  cell  system.  The  results 
of  present  study  can  be  used  in  a  comprehensive  model  of  a  PEM 
fuel  cell  to  account  for  the  effect  of  clamping  force  on  fuel  cell 
performance. 
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